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Abstract. We present bright galaxy number counts in the 
blue (16 < Bj < 21) and red (15 < R < 19.5) passbands, 
performed over 145 sq. degrees, both in the northern and 
southern galactic hemispheres. The work was conducted 
on Schmidt plates digitized with the MAMA machine, 
and individually calibrated using an adequate number of 
CCD sequences. We find a relatively large density of bright 
galaxies implying a "high" normalization of the local lu- 
minosity function. Our counts and colour distributions ex- 
hibit no large departure from what standard no-evolution 
models predict to magnitudes Bj < 21, removing the need 
for evolution of the non-dwarf galaxy population in the 
optical, out to z « 0.2. This result disagrees with that of 
Maddox et al. (1990) on the APM catalog. We show that 
the APM and similar catalogs may be affected by a sys- 
tematic magnitude scale error which would explain this 
discrepancy. 

Key words: galaxies: evolution - galaxies: fundamental 
parameters - galaxies: luminosity function - galaxies: pho- 
tometry - galaxies: statistics 



1. Introduction 



Since Hubble ( 1926| ), galaxy number counts have been 
widely used as a statistical tool for probing the distant uni- 
verse, with the hope of constraining both its geometry and 
the evolution of its content. Modern, high efficiency instru- 
ments and detectors reach impressive magnitude limits in 
the optical: Bj « 28 (Tyson |1988|, Metcalfe et al. |l995a, 



Send offprint requests to: E. Bertin (bertin@iap.fr) 
* Based in part on observations made at tiie CERGA Schmidt 
Telescope (Observatoire de la cote d'Azur, France), at Observa- 
toire de Haute- Provence (France) and at the European South- 
ern Observatory (La Silla, Chile). ESO red survey, SERC blue 
survey and CERGA plates were digitized with MAMA (Ma- 
chine Automatique a Mesurer pour I'Astronomie) which is de- 
veloped and operated by INSU/CNRS. 



Small et al. 1995] ), although on a limited area. Galaxy 
counts done at brighter magnitudes (Bj < 21) are equally 
useful; their interpretation is much less model-dependant 
because of the smaller lookback-times and weaker cosmo- 
logical effects. They thus constitute the link between mod- 
els and deeper counts by providing a normalisation of both 
space densities and colours at low redshift. 

The problem with bright galaxy counts is that they 
obviously require substantial solid angles to be surveyed 
in order to provide statistically significant samples. Un- 
til very recently, such areas could only be surveyed in a 
reasonable time using photographic Schmidt plates. Since 
the mid-eighties, fast microdensitometers coupled with 
image analysis computer programs have been employed 
to produce automatically highly complete catalogs from 
large areas of the high galactic l atitud e sky like the COS- 
MOS (Heydon-Dumbleton et al. |l989D, th e MRSP (Seitter 



1988), or the APM (Maddox et al. 1990b) galaxy surveys 



The main difficulty when dealing with photographic 
material concerns flux measurements. Many scientific is- 
sues typically require the systematic errors to be kept 
< 0.1 mag.. This is quite a difficult task to achieve on 
large scales, and can be reached only if a large number of 
galaxy standards per Schmidt plate, over the whole mag- 
nitude range of the counts, are observed to calibrate the 
data (Metcalfe et al. 1995b). Having such a high density 



of photometric standards is practically not possible with 
very large Schmidt plate surveys, because of the huge ob- 
serving time it would require. The catalogs extracted from 
these surveys are therefore likely to be hampered by sys- 
tematic errors in their photometry, rather than by limited 
statistics. 

In this paper we re-examine galaxy number counts in 
the blue and red photographic passbands over the magni- 
tude range 16 < Bj < 21 and 15 < Rf < 19.5. Our sur- 
veyed area is modest (140 sq. deg.), but has been carefully 
calibrated using a fair density of CCD standards (w 700, 
half of them beeing galaxies), as an attempt to keep pho- 
tometric systematic errors ^0.1 mag. over the whole mag- 
nitude domain of the survey. The procedure is described 
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in detail in i^, as well as the general data processing. The 
galaxy number counts are derived and compared to previ- 
ous studies in Model predictions about number counts 
and galaxy colour distributions are tested in §||. We finally 
discuss the implications of our results in m. 



2. The data 

Our optical catalog was originally built as part of a pro- 
gram of identification of faint IRAS sources on Schmidt 
plates (Bertin et al. , in preparation) near both ecliptic 
poles. The 7 selected fields are characterized by a very 
low infrared cirrus emission, which ensures that the inter- 
stellar extinction is negligible in the visible. Their coordi- 
nates are given in Table 0. To each field corresponds one 
plate in the blue photographic passband, and one in the 
red. Northern plates were taken in France at the CERGA 
0.9m Schmidt telescope, while southern ones are copies 
from the ESO and SERC surveys. 



3. Data processing 
3.1. Digitization 

The 14 plates were scanned in density mode, with the 
MAMA microdensitometer, using a 10 /xm step (0.65"). 
The MAMA is able to scan a full plate in a few hours, and 
provides images with a usable dynamic range quite large 
compared to similar machines like COSMOS or APM: 
more than 3 in density (a description of the microden- 
sitometer and its performances can be found in Berger 
et al. 1991). Because of limitations in acquisition and stor- 
age capabilities with the MAMA at that time, the image 
data were acquired as 12'xl2' frames, and recombined as 
18'xl8' images, with an overlap of 6' (which therefore de- 
fines the maximum size allowed for objects to be reliably 
measured). 

To insure a maximum reliability of the catalog, only 
the central "clean" 4.5°x4.5°of each field were kept for 
analysis (ESO and CERGA plates cover only slightly more 
than 5°x5° in total, and full coverage was not needed 
here). The exact areas retained for the final catalogs are 
reported in Table |l|. 



3.2. Astrometry 

The standard MAMA procedure was adopted to cal ibrate 
the plates astrometrically (see Berger et al. 199l|) . Be- 
tween 150 and 300 suitable astrometric standards (PPM 
catalog, Roser & Bastian 1991 ) can be found per plate. A 
third order fit on projected coordinates leads to a mean 
residual of « 0.2", similar to the value obtained when 
comparing bright object positions between the two pass- 
bands. However, this figure degrades noticeably (~1") at 
the extreme borders of some plates. We had to adjust ESO 
plates 200 and 249 to their blue counterparts SERC 200 



and 249 using a neural network mappingtl on bright stars 
to secure blue/red object matching. 

3.3. Detection 

A dedicated version of the SExtractor software (Bertin & 
Arnouts 1996, hereafter BA96) was used in PHOTO mode 
for source extraction. The standard 2 pixels (1.3") FWHM 
Sextractor's convolution mask was applied in order to im- 
prove the detectability of low surface brightnesseaS. Al- 
though emulsion noise becomes strongly correlated on 
small scales (especially on copies), this technique enables 
one to lower the detection threshold to a secure 1.5(t of the 
sky background fiuctuations, corresponding to mean sur- 
face brightnesses /ig^ = 24.9 mag.arcsec"^ (CERGA blue 
plates), ^yIy ~ ^^-^ mag.arcsec"^ (CERGA red plates), 
/igj = 25.9 mag.arcsec"^ (SERC blue plates) and //j^^ = 
23.8 mag.arcsec"^ (ESO red plates). 

The deblending capabilities of SExtractor give no limi- 
tation to the size of objects which can be extracted; spiral 
arms and other peripheral substructures of bright galaxies 
are not "separated" from the central region. Therefore the 
catalog is expected to be complete for all objects smaller 
than the overlap between subframes: 6' (corresponding 
here to Bj w 12). 



3.^. Star/galaxy separation 

Star-galaxy separation was performed with a dedicated 



neural network (Bertin 1994) trained on a sample made 
of bright objects classified by eye, and fainter ones from 
the CCD photometric fields (i |3.5.3| ). The latter were clas- 
sified automatically using SExtractor's tunable neural net- 
work classifier (see BA96). Briefly, each detected object is 
translated into a pattern vector containing basic informa- 
tion about its profile (7 isophotal areas plus the peak den- 
sity). This combination of parameters, added to the ability 
of neural networks to deal with complex distributions in 
multidimensional space, leads to a very robust star/galaxy 
classifier. In particular, objects affected by optical distor- 
tions (in the corner of the plates) or merging (we chose 
to classify undeblended pairs as "stars" if and only if the 
brightest component is a star) are properly handled. Un- 
fortunately, as we will see, such a set of parameters is not 
well fitted to the classification of very bright stars and 
galaxies, which show similar profiles. 

A minimum of 500 sample patterns is necessary for 
the neural network to achieve performances close to its 
asymptotic values (Bertin 1996); we used « 800 per field. 



^ some other examples of using neural networks as an in- 
terpolation tool in mult idimensio nal an alysis can be found in 
Serra-Ricart et al. 1995| or Bertin 



199e 



^ In SExtractor, convolution is only applied to the template 
frame on which objects are detected, not to the image itself; 
subsequent measurements are unaffected by this operation. 
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Table 1. Schmidt fields 



Field 


a (2000.0) 


5 (2000.0) 


Area* 


Plates 




h 


m 


O 


1 


sq. degrees 


"blue" 


"red" 


3'?-45° 


03 


46 


-44 


40 


21.57 


SERC 249J 


ESQ 249F 


3':'-50° 


03 


32 


-49 


40 


21.89 


SERC 200J 


ESQ 200F 


4^-4:5° 


04 


13 


-44 


40 


20.74 


SERC 250J 


ESQ 250F 


4':'-50° 


04 


01 


-49 


40 


20.38 


SERC 201J 


ESQ 201F 


9^'+78° 


09 


22 


+78 


00 


20.38 


CERGA 2970J 


CERGA 3048F 


15'? +50° 


15 


36 


+50 


00 


20.22 


CERGA 3074J 


CERGA 3063F 


16^+42° 


16 


31 


+42 


17 


20.29 


CERGA 3078J 


CERGA 3084F 



* used in the final catalog. 



Extensive visual checks and colour histograms show 
that stellar contamination among catalogued galaxies due 
to misclassifications is less than 4% for 14 < Bj < 19.5, 
rising to 10% at Bj « 20 (CERGA plates) and Bj « 21 
(UKST plates). Bj = 20 and Bj = 21 therefore define the 
upper magnitude limits for the northern and southern sets 
of plates, respectively. The estimated loss of galaxies due 
to misclassifications is similar to stellar contamination at 
faint magnitudes. We can rely on statistics from the iden- 
tifications of optically bright IRAS galaxies (Bertin et al. 



1996), which contain less than 10% of objects classified as 
"stars" for Bj > 16. As a high fraction of IRAS galax- 
ies turn out to be rather compact and difficult to sepa- 
rate from point- sources on photographic images (Suther- 
land et al. 1991 ), this gives us an upper limit of the loss 
we might expect at these magnitudes. It is more difficult 
to estimate it precisely at the bright end of our counts 
((Bj < 16), because stars outnumber galaxies by a fac- 
tor > 50. We examined 600 detections with Bj « 15.5, 
classified as stars, and found indeed that 4 of them were 
galaxies. There seems therefore to be some loss of galaxies 
at the bright end of our catalog, and counts below Bj = 16 
(or Ri? — 14.5) should only be considered as lower limits. 



case we rely on corrected isophotal magnitude. Such a sit- 
uation occurs in less than 20% of the cases in our images; 
SExtractor's total magnitude is thus essentially an aper- 
ture magnitude. 

The behaviour of SExtractor's magnitudes on simu- 
lated Schmidt plate images is shown in Fig. |l|. The pre- 
dicted fraction of flux measured for galaxies is remarkably 
constant with magnitude, even beyond the completeness 
limit. As can be seen, the mean off'set one needs to apply 
to SExtractor's magnitude to get an estimate of "total 
magnitude" is —0.06 mag. 
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3.5. Photometry 

3.5.1. Estimation of magnitudes 

The way SExtractor estimates "total" magnitudes on 
CCD images is described in details in BA96. In the 
"PHOTO" mode, the p rogram applies a density-to-intensity 
transformation (see §3.5.2) before summing pixel fluxes. 
Briefly, for each object, two kinds of magnitudes are com- 
puted. The first one is an improvement of Kron's ( 1980 ) 
method: it measures the fiux integrated in an elliptical 
aperture whose size and shape are function of the object's 
profile. The second one is isophotal, and corrects for the 
fraction of flux lost in the wings (assuming a gaussian pro- 
file) as in the APM survey (Maddox et al. |l990c| ). As this 
second method is more subject to biases than the first one, 
we take the aperture magnitude as an estimate of "total" 
magnitude, unless it is suspected to be contaminated by 
the presence of neighbours by more than 0.1 mag; in which 



Tme total magnitude 

Fig. 1. Prediction of the mean flux lost (expressed in magni- 
tude) by SExtractor's isophotal, corrected isophotal and adap- 
tive aperture magnitudes, as a function of the "true" total 
magnitude. These estimates are from simulated Schmidt plates 
images, generated as described in BA96. The simulations are 
done in the blue photographic passband; they assume a loga- 
rithmic response of the emulsion and image parameters typical 
of CERGA Schmidt plates, with a seeing FWHM of 2". The 
apparent drop seen at Bj = 16 is an artefact due to poor 
statistics. 



3.5.2. Linearization of the photographic response 

In PHOTO mode, SExtractor transforms densities D to in- 
tensities / during image analysis assuming a logarithmic 
response of the emulsion in the "linear part" : 



(1) 
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where 7 defines the contrast of the emulsion (as measured 
with MAMA), and Jo is the intensity zero-point. Both 
parameters were adjusted for each plate by minimizing 
the on "total magnitudes" estimated by SExtractor 
for a set of CCD standard galaxies. For all plates but 
one we found 2 < 7 < 3. Differential desensitization and 
telescope vignetting were compensated for by substracting 
the local background density from D before applying the 
calibration law (Eq. |l|). This correction is perfectly valid 
here as no interstellar emission is expected in these fields, 
and reduces large-scale sensitivity variations to ^ 0.1 mag 
(Bertin |l99|). 

Finally, for galaxies, a linear least-square fit was ap- 
plied separately to the magnitude scale in each field to 
correct for residual differences of ~ ±3% in the slope of 
the relation between photographic and CCD magnitudes. 
Bright stellar images are affected by a mix of photographic 
saturation, complex chemical proximity effects and, above 
all, light diffusion inside the emulsion during the scanning 
process. Stars thus need a specific calibration which was 
carried out using polynomial and neural network fits be- 
tween their photographic and CCD magnitudes. 

Eq. (0) provides a good fit to photographic calibration 
curves of sky-limited exposures measured with MAMA 
(Moreau |l992 , Bertin 1996 ), although it does not take 
into account possible saturation effects. Thanks to the 
fairly large dynamic range of MAMA, saturation affects 
our galaxy magnitudes to a lesser extent than e.g. APM or 
COSMOS scans. Nevertheless, southern copy plate images 
unambiguously show some signs of saturation on cores of 
galaxies with Bj ^ 17. We therefore expect that both sys- 
tematic and random errors increase significantly towards 
brighter magnitudes on these plates (see Metcalfe et al. 



1995b for a discussion about the consequences of satura- 



tion on galaxy photometry). 

3.5.3. CCD calibration and colour equations 

A total of 30 photometric CCD fields (w 4 per Schmidt 
field, centered on groups of bright galaxies) were taken es- 
sentially at the OHP 1.2m and the ESO 2.2 and 3.6 meter 
telescopes, and matched to Cousins B and R magnitudes. 
SExtractor was also used at all stages of the CCD data 
reduction, leading to "total magnitudes" for about 1000 
objects. As northern and southern plate sets have been 
taken with slightly different filter combinations, we natu- 
rally expect differences in their related colour equations. 
However we need to keep the survey as close as possible to 
a uniform (and as standard as possible) passband system. 
Following Blair & Gilmore ( |1982| ), we define here the blue 
passband Bj as 



Bj = B-0.28(B- V) wB-0.19(B-R) 



(2) 



which is the transformation also adopted in the APM sur- 
vey. The R_F band is defined as 



B, V and R are in the Cousins systemi. Figure ||shows the 
distribution of the residuals as a function of the Bj — R 
CCD colour index for each of the 4 types of plates. Al- 
though they are moderate, one can see systematic colour 
effects. The "best fit" colour coefficients are reported in 
Table B.5.3. For CERGA blue plates the results are in 
good agreement with the estimates of Majewski (1992), 
who finds Bi„,.,+GG385 = B - 0.23(B - V). We find a 
stronger colour coefficient for the SERC passband than 
the one traditionally used at the APM (Bjjj^^j^QQsgs — 
B — 0.28(B — V)), and confiicting even rapre with the COS 



MOS (B„:,„ 



B - 0.23(B - \)p. However our es- 



timate is in remarkable agreement with the recent and 
accurate measurements by Metcalfe et al. ( 19951:1 ), which 
give Bi„,.,+GG395 = B„,o - 0.35 (B - V) and B„,o ~ B. 

Lastly, we find in average the red photographic pass- 
bands slightly bluer than the Cousins R. This contradicts 
the traditional correction established — photoelectrically, 
and then extrapolated to photography — by Couch & 
Newell ( |198Cl| ), who proposed a response for the combi- 
nation IIIaF+RG630 redder than Cousins R. Note that 
more recent determinations (Cunow & Wargau 1993D also 
support the idea that the Couch & Newell equation is in- 
appropriate. 

In the end, the photographic magnitudes from each 
plate where corrected using the equations of Table 3.5.2 
to yield magnitudes in a unified system. Given the small- 
ness of the colour coefficients, the resulting degradation in 
photometric accuracy is negligible. The rms residual of the 
calibration above the completeness limit ranges between 
0.07 and 0.17 mag (including the contribution from large 
scale inhomogeneities). Assuming that the intensity scale 
is perfectly linearized, the formal uncertainty on the indi- 
vidual plate zero-points would be ~ 0.02 mag. But such 
an assumption is by far too optimistic with photographic 
plates. Given the number of bright standard galaxies per 
plate (between 1 and 3 per magnitude), and the rms uncer- 
tainty on their magnitudes (about 0.15 mag), we estimate 
the individual zero-point systematic errors to be ^ 0.1 
mag in the range 15 < Bj < 21, and 14 < Rj- < 19.5. As 
Figure ^ shows, below these fluxes the photometric errors 
grow rapidly. 

3.5.4. Comparison with other photometries 

Figure || shows the differences between photo and CCD 
magnitudes for our CCD standard galaxies, as well as 
the comparison with the photometry of bright galaxies 



Although our Rf is photometrically equivalent to Cousins 
R within the measurement errors, throughout this paper we 
denote red photographic magnitudes "Rf" instead of "R". 
^ It seems that through the years some confusion has often 
been made between the IIIaJ-|-GG385 and IIIaJ-|-GG395 com- 



R^ = R 



(3) 



binations (see for instance Shanks et al. 1984), and yet these 
prove to lead to somewhat different (although close) passband 
definitions. 
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Fig. 2. Differences between "pure" (but calibrated) photographic and CCD magnitudes as a function of the B j — R colour index 
for the brightest standards: Bj < 20 (CERGA blue), Bj < 21 (UKST blue), R < 18 (CERGA red) and R < 19 (ESO red). 
The standards contain about 2/3 of stars and 1/3 of gala^xies (which show similar behaviour). A few points (crosses) deviate by 
more than 3a from the mean relation and have not been taken into account for the chi-square fit (solid line). 



Table 2. Photographic passbands and their link to standard photometric systems 



Plate type 


Emulsion+filter 


Colour equations 




CERGA "blue" 


IllaJ + GG385 


Bj + (0.033 ±0.015) (Bj 


-R) 


B - 0.24(B - 


-V) 


UKST "blue" 


IllaJ + GG395 


Bj - (0.060 ± 0.023) (Bj 


-R) 


B-0.35(B- 


-V) 


CERGA "red" 


IIIaF + RG610 


R+ (0.044 ±0.017) (Bj - 


-R) 


R± 0.036(B 


-R) 


ESO "red" 


IIIaF + RG630 


R± (0.048 ±0.019) (Bj - 


-R) 


R± 0.039(B 


-R) 



by several authors. Nine of the galaxies measured by Met- 
calfe et al. (1995b) are found in our catalog; we find mean 
differences (MAMA-Metcalfe) ABj = -0.04 ± 0.08, and 
AKf = —0.07 ± 0.05. Note that 3 of these galaxies also 
figure among our CCD standards; in both B and R, our 
magnitudes are in agreement with theirs within 0.02 mag. 

Although their photographic photometry does not 
reach the same accuracy and is more subject to biases 
than CCD or photoelectric standards, galaxies from the 
LV catalog (Lauberts & Valentijn 1989| ) have proven to 
be in good agreement with the RC3 system (Paturel et al. 



1994), and can be used to trace an eventual large sys- 
tematic trend at bright magnitudes (Bj < 16). Over the 
4 southern fields we find mean differences (MAMA-LV) 
ABj = +0.10 ± 0.05 and ARj. = ±0.06 ± 0.05. These 



discrepancies might be interpreted as some loss of flux 
at the bright end (Bj « 15) of our catalog. However, in 
the same magnitude range, "CCD galaxies" photometered 
here or by Metcalfe et al. do not show this offset, which 
could argue for a small, local zero-point error in the LV 
magnitudes. 

Some bright galaxies (Bj ^ 15) from the northern 
plates do also have photoelectric B and V "total" mag- 
nitudes in the RC3 catalog. From these we get a mean 
difference (MAMA-RC3) ABj = ±0.15 ± 0.08 which re- 
veals the influence of plate saturation at the bright end of 
our catalog. 

Finally, 5 of the galaxies photometered with a CCD 
by Maddox et al. ( |1990cD to calibrate the APM survey 
lie in our field 201 and give an offset (MAMA-APM) 
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ABj = -0.16 ± 0.14. More than this significant offset, 
the unexpectedly large dispersion of magnitudes — almost 
0.3, that is, as much as with the photographic LV sample! 
— casts some doubt over the reliability of the Maddox 
et al. calibration set for this plate. 

In conclusion, we believe our magnitude scale to be 
free from any large systematic error, except brightwards of 
Bj « 15, where fluxes might possibly be underestimated 
by « 0.1 mag. 

3. 6. Merging of catalogs 

The "blue" and "red" catalogs were cross-identified in al- 
pha, delta, to yield a unique two-color catalog. Differences 
in seeing or image quality, as well as imbricated detec- 
tions (i.e. stars lying on disks of galaxies) were handled 
through a complex matching algorithm taking into ac- 
count the shape of detected objects. This procedure leads 
to a highly reliable catalog, virtually suppressing all false 
detections like pieces of hair, satellite trails, optical ghosts 
or spikes around bright stars (these are generally classified 
as galaxies). One might fear some loss of objects with ex- 
treme colours in the final catalog; however the fraction of 
non-paired detections is almost constant with magnitude, 
and is about 1% (most of which are spurious), rising to 
5% in the last half-magnitude bin imposed by star/galaxy 
separation. 

4. Galaxy number counts 

Results 

The differential number counts are shown for each field 
in FigJ^. One can notice a great homogeneity over the 
whole magnitude range between the plates, except for the 
16'.'-f42° field. In fact this field contains no less than 12 
Abell clusters, and is part of the Hercules supercluster. 
One can estimate the variance in number counts at a 
given depth on a Schmidt plate from the Maddox et al. 
( |1990aD angular two-point correlation function (e.g. Pee- 
bles 1980). We find that the projected density of galaxies 
in the 16l'-|-42° field exceeds the mean one from the other 
fields by a factor of 4 times the expected deviation at 
Bj = 15, reducing to « 2 for Bj > 17. This field is there- 
fore one of the densest in the sky at the bright end of 
our catalog, and including it in the total number counts 
shown in Fig.|| severely changes the slope at Bj <^ 16. 
For this reason we decided to discard it from the samples 
for the statistical analysis of number counts, remembering 
however that we might underestimate in this way the true 
density of bright galaxies. 

4-2. Comparison with previous counts 

Bright galaxy counts published so far were all conducted 
on photographic material, in quite similar passbands. This 
enables us to compare easily our results with previous 



work. The transformations we applied are described be- 
low. 

M any o f the existing brig ht ga laxy counts (Shanks 
et al. 



1984 



1986 Heydon-Dumbleton 



Stevenson et al. 
et al. |1989| ) mak e use of the bj passband definition estab- 
lished by Kron ( 1980| ), and which transforms to our Bj 
through 



Bj = bj - 0.05 



(4) 



assum ing (B — V) = 1.0. The recent counts of Weir et al. 
( 1995 ) were done in the Gunn-g passband, for which they 
give the approximate transformation 



B. 



5-0.5 



(5) 



to the APM survey magnitude system (Maddox et al. 
1990c), which is identical to ours. 

G unn-r observations (Sebok 1986| , Picard 1991b| , Weir 
at al. 1995) were transformed to the Cousins R system 
using 



R = r + 0.341 



(6) 



from a fit to 16 Landolt standards, giving a standard er- 
ror of only 0.014 mag (P.Fouque, private communication). 
Shanks et al. (1984) rp magnitudes were converted with 



R = rF -0.09 



(7) 



inserting (B — R) = 1.5 in their equation (5). 

As shown on Fig. our counts agree well with those 
from other studies, especially with the preliminary counts 
done on POSSII plates (Weir et al. 1995 ), despite the un- 
certainties in the conversion from their photometric sys- 
tem to ours. In R, the major discrepancies observed are 
with Shanks et al. ( |1984D , and the Picard ( |l991b[ ) northern 
and southern counts, although our results lie in-between. 
The Schanks et al. counts are based on only one Schmidt 
plate, so the fact that they show an offset is perhaps not 



surprising. As pointed out by Weir et al. (1995), it is yet 
unclear whether the difference between the Picard counts 
originates from the presenc e of v ery large scale struc- 
tures, as claimed by Picard ( 1991a ), or are simply caused 
by unexpected photometric errors in his photometry. We 
note here that one of our plates lies at about 20°from Pi- 
card's northern field, and that no particular enhancement 
in galaxy density is seen there. 

At their faint end (20 < Bj < 21), our b lue co unts 
are about 10% lower t han those of Weir et al. ( 1995| ) and 
Maddox et al. (1990d); however the raw counts (Fig. H) 
do not indicate any obvious drop in completeness at that 
level. As a matter of fact, two points should be consid- 
ered here: (1) half of the difference (5%) can be explained 
by detections that were not matched between the two 
colours and dropped in the final catalog (although there 
is no proof that these are real objects), and (2) Weir et al. 
(1995) showed on simulations that both their counts and 
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Fig. 4. "Raw" differential number counts in Bj and Rf for each of the Schmidt fields (continous line = southern fields; dashed 
lines — northern fields). Note the excess of galaxies in the northern field 16^'-|-42° (thicker lines), especially at bright magnitudes. 



those derived by the APM survey might be slightly over- or Heydon-Dumbleton et al. ( |1989| ). The recent Weir et al 



estimated (by about 10%) for Bj ^ 20. Therefore it is 
likely that this difference at the faint end between our 
counts and others should not be interpreted as something 
significant, but only as a consequence of a different data 
processing. This is confirmed by the recent medium-deep 
CCD counts by Arnouts et al. 1996, which are in perfect 



agreement with ours over this magnitude range, in both 
colours. 

But the most interesting discrepancy in this compar- 
ison is the one seen with some other galaxy counts at 
bright magnitudes. At Bj = 16.5 for instance, we count 



(1995) counts stay however in perfect agreement with ours 
down to their brightest limit. As the APM survey, because 
of its statistical weight, is traditionaly used to normalize 
galaxy counts models at bright magnitudes, it is worth 
investigating further what may be the origin of this dif- 
ference. Recently, Metcalfe et al. ( |1995b[ ) have discovered 
some scale error on the 17 ^ Bj < 18 domain in the APM 
galaxy magnitudes. The APM survey shares four Schmidt 
fields in common with us, and the corresponding part of 
the catalog, down to Bj = 21 was kindly provided by J. 
Loveday for comparison. 



nearly twice as many galaxies as do Maddox et al. (1990d) 
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Bj Rp 

Fig. 5. Total differential number counts in Bj and Rf in our catalog, including (open circles) or not (filled circles) the 16'?+42° 
field. Error bars indicate poissonian (Ict) uncertainties whereas lines bracket the rms excursion of the total counts assuming the 
galaxy density is uncorrelated from plate to plate. 



Figure |^ shows the difference between the APM cat- 
log and our magnitudes for each of these four fields 
Although both magnitudes are in agreement within 
w 0.1 mag. at the faint end (Bj > 20), one can imme- 
diately notice a systematic difference at brighter magni- 
tudes, reaching 0.3 to 0.5 mag. at Bj « 16—17. To find out 
whether this was due to our calibration or was intrinsic 
to the APM data, we cross-identified our list of south- 
ern standard galaxies with the APM catalog to produce 
a plot similar to Fig. |^. Figure ^ shows indeed that there 
seems to be some problem in the magnitude scale of this 
APM galaxy sample. The magnitude overestimate starts 
to rise at Bj « 19.5 and increases until Bj « 17 where it 
culminates at about 0.4 mag. Brighter than this point, it 
is unclear whether it stays at the same level or decreases 
again, but at Bj w 15 — 16 the offset is still at least about 
0.2 mag. Of course these estimates are in principle only 
valid for the small subsample of the APM catalog consid- 
ered here, and might not be applicable to the full catalog. 
But one can see that compensating for this average trend 
in the magnitude scale would bring the total APM counts 
in very good agreement with ours. 



^ The version of the APM catalog shown here uses plate- 
to-plate matching based only on the faintest magnitude bin 
(19.5 < bj < 20.5). The actual version used for statistical 
studies, limited to bj > 17, uses matching in three mag- 
nitude ranges (G. Dalton, private communication). However, 
a comparison between the two versions for our zone shows 
an average difference (unmatched— matched) peaking to only 
~ -1-0.06 mag. at bj = 17. 



5. Comparison with models 

5.1. Numbers counts and the normalization of the "local" 
luminosity function 

A recurring problem with previous counts has been the 
normalisation of models at bright magnitudes. Normal- 
izing passively evolving models at Bj = 17 leads to an 
apparent excess of already 50 to 100% at Bj « 20 — 21. 
This apparent phenomenon is hard to reconcile with the 
redshift distributions to Bj ^ 23 — 24 (Glazebrook et al. 



1994, and references therein), which suggest a higher nor- 



malisation of the field luminosity function. As our data 
"push" the bright end of galaxy counts to higher values, 
we can now test how local determinations of the luminos- 
ity function are able to match the counts over the domain 
16 < Bj < 21. 

We therefore constructed a very simple non-evolving 
model of galaxy counts, intended to be valid at least to 
Bj ^ 20. K-corrections were computed by integr ating the 
spectral energy distributions from Pence ( 1976| ) through 



the bj and rp photographic passbands given by Couch 



Newell ( |1980D , which are very close to our blue and 
red passbands. The morphological mix of galaxy types 
was taken from Shanks et al. ( 1984 ). This model does 
not include luminosity or number evolution as in deeper 



count models (e.g. Guiderdoni & Rocca-Volmerange 1990 



McLeod & Riecke 1995), but is sufficient for our normal- 



isation purpose. We chose an Einstein-de Sitter universe 
for simplicity, although at this depth {z ^ 0.2) the counts 



E. Bertin and M. Dennefeld: Galaxy evolution at low redshift? - I. Optical counts 



9 



1000 



CO 



t3 



100 



10 



0.1 



~i — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



ft 



*6 



9 

6 



yf- ft 



• This work 

□ Stevenson et al. (1986) 

X Heydon-Dumbleton et al. (1989) 

O Maddox et al. (1990) 

* Weir et al. (1995) 



1" I I I I I I i I I \ \ 1 I I I L 



1000 



tiO 
CO 



I 

T3 



14 



16 18 



20 



22 



100 



10 



0.1 



1 


' ' 1 


1 1 1 1 1 


1 1 1 1 1 1 


1 


— 






• 

* 


-= 








A 
Ok A 




- 

- 






J, s 


- 
- 


— 








— 
z 


- 








- 




















"= 
















s 








i 




• This work 










□ Shanks et al. (1984) 






i 




O Sebok (1986) 










A Picard (1991) north 






i 




A Picard (1991) south 










* Weir et al. (1995) 




~ 1 


, , 1 


1 1 1 1 1 


, 1 , , , 1 


1 ~ 



12 



14 



16 



18 



20 



R 



Fig. 6. Comparison of published bright g alax y counts converted to our Bj,R photometric system. Arrows indicate counts 
underestimated because of incompleteness (3.4). 



< 
< 

I 

< 



-1 - 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




UKST 200 
I I I I I I I I I I I I 



< 

< 

0.5 



m 

I 



12 



14 



16 18 
Bj (MAMA) 



20 




16 18 
Bj (MAMA) 



1 - 







g-0.5 

^ - 1 P 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




UKST 201 

I I I I I I I I I I I I 



22 



12 14 16 18 20 22 
Bj (mama) 



< 1 
0.5 h 



m 

I 







£-0.5 
< 



1 - 



~i — I — r 



~i — I — I — I — I — I — I — I — I — r 




UKST 250 



22 



12 



14 



16 18 
Bj (mama) 



20 



22 



Fig. 7. Comparison of galaxy magnitudes between the APM and our MAMA catalog on four UKST Schmidt plates. The tilted 
frontier on the right side of the plot is caused by the Bj — 21 limit in the APM subsample. 



10 



E. Bertin and M. Dennefeld: Galaxy evolution at low redshift? - I. Optical counts 



prove to be quite insensitive to (reasonable) cosmological 
parameters (e.g. Yoshii & Takahara 1988 ). 

Figure ^ shows our data compared to simple galaxy 
counts models in the blue photographic passband with- 
out evolution as described above, using two different lu- 
minosity functions (LF). The first one is from Efstathiou 
et al. ( |1988 , hereafter EEP), from a compilation of re- 
sults obtained on several redshift surveys, and t he second 
one, more recent, is from Loveday et al. (1992, hereafter 



LPEM) from a sparse redshift survey of APM galaxies. 
Schechter parameters are given for both in Table A 
third, colour dependant, LF from Shanks (199C), with the 
normalisation adopted by Metcalfe et al. (1991) is also 
shown for comparison. As can be seen, the EEP luminos- 
ity function fits the data much better than the one from 
LPEM, which are however, as expected, in good agree- 
ment with the APM counts for Bj < 17. In fact, most of 
the discrepancy between the two luminosity functions lies 
in the value of M*^ , with a difference close to 0.4 mag, 
i.e. about the same one gets when comparing APM mag- 
nitudes to those of standard galaxies at Bj w 17 (Fig. ||). 
The best agreement is found with the Shanks LF model. 
This is not a surprise, as the two other models do not dis- 
tinguish the luminosity functions of early and late-types, 
which are indeed quite different; and this shows up when 
K-corrections become important. 
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Fig. 8. Same as Fig. ||, but with the APM catalog. The dashed 
line marks the magnitude limit of the APM subsample. 

Figure ^indicates that the norma lisation of the Shanks 
LF adopted by Metcalfe et al. ( 1991 ) is in excellent agree- 
ment with our own counts, but the latter are still slightly 
too steep on their bright side compared to any non- 
evolving model. Can this be interpreted as evolution? One 
can hardly conclude on the basis of the data presented 
here, as a 10% misclassification of galaxies added to a 
0.1 mag. systematic error in measured fluxes at Bj w 17 
would be enough to produce this effect. 



1000 




Fig. 9. Comparison of our Bj galaxy counts (filled circles) and 
those of the APM survey (open circles) with the simple models 
described in the text. Errorbars are la uncertainties deduced 
from the field-to-field scatter within each bin (assuming galaxy 
density is uncorrelated from one plate to another). Arrows in- 
dicate points affected by incompleteness. The 3 curves corre- 
spond to different lumin osity f unctions used in the no-evolution 



model: Efstathiou et al. 



1988 



(upper solid line), Loveday et al. 



1992 (lower solid line), and Shanks 1990 as normalized by Met 



calfe et al. 1991 (dashed line) 



5.2. Galaxy colours 

Galaxy colour distributions for the total catalog (including 
the field 161^4-42°) are presented in Fig. |o[ They are in 
good agreement with the compilation presented by Koo 
& Kron ( |1992| ), if one makes allowance for the passband 
differences between our system and theirs. 

We reemployed the Shanks LF model with no evolution 
described above to model the observed colour distribution. 
Rest-frame galaxy c olour s per Hubble type were taken 
from Metcalfe et al. (1991) and converted from their CCD 
photometric passbands to our photographic system with 
(Bj - Rf) = 0.91(6 - R)ccD. Ingredients of the LF are 



gathered in Table 5.2. Photometric errors were included in 
the model, adopting an error distribution with magnitude 
determined in the overlap between catalogs (see Bertin 
1996). Photometric errors do not only smear out details 



in the observed colour distribution, they also distort the 
wings of the distribution for the dimmest objects (that is, 
when the error strongly evolves with colour). As our red 
plates are here significantly less sensitive than the blue 
ones, the two lower graphs in Fig. |l^ exhibit a shallower 
tail on their blue side than on their red side. 
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Table 3. Global Schechter parameters of local luminosity functions (_ffo ~ 100km. s ^Mpc ^) 





Q 


A/*, 


0* (Mpc^^) 


Efstathiou et al (ll98j) 
Loveday et al. (|l992D 


-1.07 ±0.05 
-0.97 ±0.15 


-19.97" ± 0.10 
-19.50 ±0.13 


(1.56 ±0.34) X 10~^ 
(1.40 ±0.17) X 10-'^ 



Converted from Bt using Eq. (4. Id) from Efstathiou et al. 



Observed and predicted distributions prove to be in 
very good agreement, except in the brightest subsample, 
where we expect some problems with the saturation and 
the robustness of star/galaxy separation. One can notice 
a small (« 0.1 mag) systematic colour offset between the 
two, which might be imputed to uncertainties in our pho- 
tometric modelling. But no obvious evolution of colour 
with magnitude is detected, within the uncertainties, up 
to Bj < 21. 



Table 4. Luminosity function parameters for the colour dis- 
tribution model {Ho = 100km.s"^Mpc"^) 





a 


Mb, 


(j)* (Mpc" 




(Bj - Rf) 


E/SO 


-0.7 


-19.6 


1.2 X 10" 


-3 


1.32 


Sab 


-0.7 


-19.6 


3.9 X 10" 


-3 


1.15 


Sbc 


-1.1 


-19.9 


6.9 X 10" 


-3 


0.96 


Scd 


-1.5 


-20.1 


2.7 X 10" 


-3 


0.70 


Sdm 


-1.5 


-20.1 


1.7 X 10" 


-3 


0.62 



for B j <^ 17, or over the full sky for Bj ^ 12 (Paturel 
et al. 1994 ), the an alysis of a large spectroscopic subsam- 



ple (Loveday et al. 1992 ) seems to exclude this possiblity 
in the APM survey. 

Statistics of CCD calibrations (Maddox et al 



1990d, 



Loveday et al. 1992) done on the APM magnitudes do not 
reveal any important trend like the one in Fig. ^ and ^ 
This is quite puzzling, and one might wonder if the effects 
seen here are not purely local and compensated somewhere 
else in the APM survey. If this were true, correlated sys- 
tematic magnitude errors of ~ 0.4 mag., over such large 
areas (100 sq. deg. in our case) would artificially boost the 
two-point correlation function on scales <; 5° in Bj ^ 18 
magnitude bins. Th is does not appear prominently in Fig. 
2 of Maddox et al. (1990a). Another clue is that the trend 



foimd by Metcalfe et al. ( 1995a ) concerns the whole area 
of the APM catalog, and therefore supports the hypothesis 
of a global magnitude scale error (which does not exclude 
large scale variations of the zero-point). 



6. Discussion 

6.1. A magnitude scale error in the APM data? 

The apparent rapid increase in number counts of galaxies 
in the magnitude range 16 ^ Bj ^ 19 reported by Mad- 
dox et al. ( 1990cj) i s thus absent from our counts (or those 
of Weir et al. 1995), although our slope of dN/dm is still 
slightly higher than what predict no-evolution models at 
Bj < 17. As the comparison with the APM data seems 
to indicate, this might be essentially due to a difference 
in magnitude scale of « 0.1-0.2 mag. /magnitude. Serious 
suspicions about systematic residual errors in the APM 
magnitudes at Bj » 17 — 18 have recently been raised 
by Metcalfe et al. ( 1995a ). If our data are correct, they 



support and even strengthen this hypothesis. Other argu- 
ments also consolidate this idea. 

First, our data (Fig. ^0|), as well as previous ones (see 
Koo & Kron 1992| ), do not present evidence for evolution 
of galaxy colours over the range 17 < Bj < 21. and rule 
out strong luminosity evolution scenarios in the visible 
(McLeod & Rieke 1995) for bright field galaxies, out to 
z « 0.2. 

One might also think of a "giant local void" that would 
affect the bright end of the counts. Although there are 
indications of large fluctuations in projected galaxy den- 
sity at the level of a Schmidt field (Fig. |, Picard |l991bD 



6.2. Uncertainties in the "local" field luminosity function 

To partially bypass uncertainties in the optimum 
Schechter parameters of the local luminosity function for 
field galaxies, count models are generally "renormalized" 
by adjusting the Schechter density parameter (jj* to counts 
at some bright magnitude. However, the Al* parameter 
itself is directly affected by any systematic error in the 
magnitude zero-point of a redshift survey. Such errors are 
known to exist at bright magnitudes in photographic cat- 
alogs (Metcalfe et al. 1995b , Yasuda et al. 1995 ) on which 
are essentially based all the determinations of the local 
field luminosity function. In fact, decreasing the M* found 
by LPEM by about 0.4 mag., like our photometric com- 
parison with the APM catalog suggests, brings their lu- 
minosity function in good agreement with both the one 
from EEP and our bright galaxy counts (Fig. |^), with- 
out having to increase (f>* . One should also note that such 
a "bright" normalisation in M* of the luminosity function 
removes any significant luminosity evolutioncl of luminous 
galaxies in the blue band to z « 0.2 (Lonsdale & Chokshi 
1991. 



® The related case of evolution in the far-infrared as seen by 
IRAS will be addressed in a forthcoming paper (Bertin et al., 
in preparation). 
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Fig. 10. Colour histograms of galaxies compared with the no-evolution model (curves) described in the text, in 4 magnitude 
intervals. The colour distribution of stars (thin histogram) is plotted for comparison in the faintest subsample. 



7. Summary and Conclusions 

We have presented a new photometric survey of bright 
galaxies with 15 < Bj < 21 and 14 < Rj^ < 19.5. Com- 
parison with our own CCD galaxy standards, as well as 
others from the literature, allowed us to estimate system- 
atic errors in the photometry to be ^ 0.1 mag over the 
whole magnitude range of the survey. Our data reveal no 
evidence of strong galaxy evolution as had been reported 
by Maddox et al. (|l990dD with the APM survey, although 
the slope of our number counts is still somewhat higher 
than what is predicted by no-evolution models brightward 
Bj « 17. A comparison with the APM data and other 
arguments suggest that the APM catalog is affected by 
a large magnitude scale error, underestimating by about 
0.4 mag the flux of galaxies with Bj ^ 17. This would 
then justify the choice of a "high" normalization of the 
field luminosity function (concerning M* as well as (j)*) 
already used by several authors to improve the fit of mod- 
els at Bj > 19 (e.g. Metcalfe et al. |l99iyi9954 Glazebrook 



coming digital large scale surveys like the SDSS (e.g. Kent 



1994) should provide this possibility. 



et al. |1994| ). 

But, given the difficulties inherent to photographic 
photometry, it would be good to confirm our analysis with 
linear detectors in order to conclude definitively. Forth- 
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